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STATEMENT  OF  THE  PROBLEM 


The  response  of  a  long  rod  penetrator  striking  armor  is  of 
considerable  military  interest.  The  complete  large-deformation  three- 
dimensional  response,  including  fracture  and  erosion,  is  difficult  and 
expensive  to  model.  The  problem  addressed  here  was  to  determine  if  the 
part  of  the  projectile  behind  the  tip  could  be  modeled  in  an  approximate 
way  using  beam  theory. 


SUMMARY  OF  KEY  RESULTS 

The  two  most  significant  results  of  this  research  are: 

(1)  Development  of  a  theory  for  treating  the  thrust- 
dominated  plastic  thrust moment  Interaction  in  a  beam 
during  loading  and  partial  unloading. 

(2)  Insight  to  the  local  loading,  unloading,  and  reloading 
plastic  response  of  a  beam  loaded  axially  and  trans¬ 
versely  at  the  end. 

A  theory  was  developed  for  modeling  the  dynamic  elastic-plastic 
response  of  a  beam-column  under  combined  axial  and  lateral  loading.  The 
theory  is  based  on  the  Timoshenko  beam  equations  and  a  new  formulation 
for  the  constitutive  relations  in  the  plastic  regime.  The  constitutive 
relations  rely  on  the  simplifying  assumption  that  the  state  of  stress  is 
dominated  by  the  thrust  as  is  the  case  for  long  rod  penetrators  impact¬ 
ing  at  small  angles  of  attack.  This  simplifying  assumption  has  three 
implications:  (1)  bending,  both  in  loading  and  unloading  at  a  local 
cross-section,  is  treated  as  a  perturbation  on  the  thrust;  (2)  the 
degree  of  work  hardening  at  a  given  cross  section  depends  only  on  the 
thrust  at  that  cross  section;  and  (3)  the  plastic  shear  strain  can  be 
neglected. 

The  analysis  is  based  on  the  Timoshenko  beam  equations.  For 
elastic  response,  the  thrust -longitudinal  strain  and  moment -curvature 
relationships  are  linear  and  uncoupled.  For  plastic  response  during 
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loading  and  unloading  over  a  portion  of  the  cross  section  (partial 
unloading),  the  resultant  thrust  and  moment  at  that  cross  section  are 
coupled.  These  resultants  are  found  by  applying  the  appropriate  stress- 
strain  relation  (elastic  or  plastic)  at  each  point  of  the  cross  section 
and  Integrating  over  the  cross-sectional  area.  A  single  parameter  can 
be  used  to  specify  the  location  of  the  boundary  between  the  elastic  and 
plastic  regions.  This  approach  was  used  to  formulate  the  stress-strain 
relations  for  loading  and  partial  unloading,  both  of  which  are  dominated 
by  thrust.  The  equations  were  solved  by  the  method  of  characteristics. 

We  attempted  to  calculate  the  response  of  a  penetrator-type  beam 
under  a  predominantly  axial  end  load  such  that  the  dominant  loading  due 
to  thrust  would  be  accompanied  by  partial  unloading  due  to  bending.  For 
small  values  of  the  angle  between  the  direction  of  the  resultant  load 
and  the  beam  axis  only  loading,  predominantly  axial  with  small  bending, 
was  produced.  As  this  angle  increased  to  values  as  small  as  8  degrees, 
partial  unloading  caused  by  the  Increased  bending  was  produced  over  a 
portion  of  the  length  of  the  rod.  However,  this  portion  soon  became 
reloaded . 

We  were  unable  to  generate  partial  unloading  at  any  location 
without  significant  reloading.  We  attribute  this  effect  to  the  fact 
that  partial  unloading  occurs  along  a  steeper  slope  of  the  stress-strain 
curve  than  the  slope  for  plastic  loading.  Therefore,  as  unloading  takes 
place  there  is  an  Increase  In  the  effective  local  beam  stiffness  resist¬ 
ing  the  superposed  bending.  This  disturbance  propagates  away  from  the 
Initial  point  of  partial  unloading  and  causes  reloading  at  nearby 
points. 

We  conclude  that  any  model  of  penetrator  response  must  include  at 
least  one  cycle  of  loading,  partial  unloading  and,  reloading. 
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ABSTRACT 


A  one -dimensional  structural  mechanics  model  Is  presented  for 
analysing  the  elastic-plastic  response  of  rods  to  slightly  Inclined  end 
thrusts.  The  theory  Is  kept  simple  by  restricting  the  development  to 
cases  In  which  the  axial  thrust  dominates,  the  bending  moment  being 
considered  as  a  perturbation  that  may  grow.  Partial  unloading  of  a 
cross  section  caused  by  excessive  bending  is  Included  in  the  treatment 
but  requires  further  numerical  study.  Solutions  are  obtained  using  the 
method  of  characteristics.  The  computer  code  describing  the  elastic- 
plastic  wave  propagation  Is  simple  and  may  be  useful  for  preliminary 
design  of  rod  penetrators. 
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1.  INTRODUCTION 


Long-rod  penetrators  for  defeating  armor  have  been  extensively 
analyzed  by  finite-difference  computer  codes  to  obtain  their  axi- 
symmetrlc  elastic-plastic  response  to  Impact.  If  the  end  loads  result 
in  stress  waves  that  are  not  axisymmetric,  the  continuum  code  solutions 
become  three-dimensional  and  complex.  As  an  aid  to  preliminary  analyses 
of  rods,  we  present  a  one-dimensional  structural  mechanics  approach  that 
is  suitable  for  treating  elastic-plastic  wave  propagation  along  rods  for 
cases  involving  small  departures  from  axlsymmetry. 

In  this  simple  theory,  we  assume  that  plane  cross  sections  of  the 
rod  always  remain  plane,  which  allows  us  to  formulate  the  problem  in 
terms  of  the  axial  and  transverse  velocities  of  the  center  of  the  cir¬ 
cular  cross  section,  the  angular  velocity  of  the  planar  cross  section 
about  the  appropriate  diameter,  and  the  stres8  resultants  consisting  of 
the  thrust,  bending  moment,  and  shear  force  acting  on  the  section.  The 
derivation  of  the  governing  equations  is  based  on  the  Timoshenko  theory 
of  beams  [1).  This  well-known  theory  takes  into  account  rotary  inertia 
and  overall  shear  distortion  of  a  beam  element.  However,  our  main 
reason  for  using  this  theory  as  a  basis  is  that  the  resulting  equations 
are  hyperbolic  and  lend  themselves  to  solution  by  the  method  of 
characteristics  [2]. 

Colton  [3]  has  used  this  structural  mechanics  approach  to  analyze 
the  elastic  response  of  earth  penetrators,  which  may  be  considered  long 
rods.  Good  agreement  was  found  between  predicted  and  experimental 
strains.  Ranganath  and  Clifton  [A]  have  also  used  this  approach  suc¬ 
cessfully  to  analyse  the  elastic-plastic  response  of  an  infinite  beam  to 
transverse  impact . 

He  have  intentionally  restricted  the  theory  so  that  the  axial 
thrust  is  supposed  to  dominate  the  bending  moment.  The  axially  applied 
end  load  monotonically  Increases  and  it  Is  kept  at  a  small  angle  of 
inclination  to  the  rod  axis.  It  is  also  assumed  that  the  shear  force 


does  not  participate  in  the  yield  condition*  These  restrictions  allow 
us  to  assess  the  usefulness  of  the  elastic-plastic  theory  in  its  sim¬ 
plest  form,  after  which  the  effects  of  relaxing  restrictions  can  be 
studied* 

Figure  1  illustrates  the  effect  of  a  dominant  thrust  with  the 
bending  moment  regarded  as  a  perturbation  that  can  grow.  Figure  la 
shows  the  end  thrust  NQ  inclined  at  a  small  angle  a  to  the  rod  axis. 

Let  ACB  be  a  cross  section  of  the  rod  where  A  and  B  are  the  top  and 
bottom  fibers  and  C  is  the  fiber  on  the  horizontal  diameter  or  mid¬ 
plane.  If  the  curvature  at  the  section  has  its  concave  side  on  top,  the 
compressive  strain  at  A  is  larger  than  the  compressive  strain  at  C  while 
the  strain  at  B  is  smaller  than  the  strain  at  C.  Also,  while  loading 
increases,  all  strains  Increase  provided  the  thrust  dominates.  The 
stress-strain  states  at  A,  C,  and  B  are  shown  in  Fig.  lb,  with  the 
attached  arrows  indicating  that  the  strains  are  increasing.  At  B,  the 
bending  that  is  causing  the  stress-strain  state  to  lag  behind  may  grow 
enough  to  cause  strain-rate  reversal,  in  which  case  partial  unloading  of 
the  cross  section  occurs.  This  eventuality  is  in  fact  Included  in  the 
theoretical  development  of  the  next  section,  in  an  attempt  to  prepare 
for  extension  of  the  simplest  form  of  the  elastic-plastic  theory.  It  is 
also  Included  to  show  how  partial  unloading  Introduces  new  wave  veloc¬ 
ities,  as  in  the  work  of  Ranganath  and  Clifton  [4]. 

2.  GOVERNING  EQUATIONS 

Figure  2a  shows  a  rod  element  with  its  attendant  axial  force,  shear 
force,  and  bending  moment,  which  at  an  axial  position  x  are  represented 
by  N,  Q,  and  M.  Figure  2b  shows  the  axial  and  transverse  displacements, 
€  and  t),  the  rotation  <t>,  and  element  overall  shear  strain  y,  all  at 
location  x.  Thus,  the  three  equations  of  motion  may  be  written  as 


MI  +  Q  -»!♦« 


(2.U) 


(2.1b) 
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Qx  -  p  a  ntt 


(2.1c) 


2 

where  p  is  the  density  of  the  rod  material,  A  -  xa  Is  the  cross* 
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sectional  area,  and  I  »  xa  /4  is  the  second  moment  of  area.  Subscripts 
x  and  t  denote  partial  differentiation. 

Because  we  assume  that  the  shear  strains  remain  elastic,  we  may 

take 

Q  -  k*  A  G  y  (2.2a) 

where  kf  Is  the  cross  section  shear  coefficient  [5]  and  G  is  the  modulus 
of  rigidity.  For  a  circular  cross  section,  we  take  kf  -  0.827  [6J.  In 
terms  of  the  rod  element  slope  and  rotation,  the  section  shear  strain  Is 

Y  -  tix  -  ^  (2.3a) 

When  the  compressive  strains  are  Increasing  over  the  entire  cross 
section,  as  shown  in  Fig.  lb,  the  Incremental  forms  for  the  generalized 
stress-strain  relationships  are 

N  -  H  A  et  (2.4a) 

M  ■  H  I  ic  (2.4b) 

where  e  -  £  Is  the  cross  section  midsurface  strain  and  ic  -  <b  is  the 
x  x 

curvature  of  the  rod  element.  The  modulus  H(e)  is  a  constant  only  when 
the  strain  Is  less  than  the  yield  strain,  the  case  when  H  becomes 
Young's  modulus  E. 

In  terms  of  the  axial,  transverse,  and  angular  velocities  u,  v, 
and  a)  given  by 

u  -  Zt  v  -  i)t  w  -  *t  (2.5a) 
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equations  (2.1a),  (2.1b),  (2.1c),  (2.2a),  (2.4a),  and  (2.4b)  can  be  put 
In  the  form 


p  A  u  -  N  *0 

t  x 

N  -  H  A  u  -  0 

t  x 

old)  -  M  «  Q 

r  t  x 

M  -Hid)  -  0 

t  x 

p  A  vt  -  Qx  -  0 

Qt  -  k'  A  G  •  -k'  A  G  d)  (2.6a) 

To  allow  the  theory  to  treat  the  case  of  partial  unloading,  which 
starts  with  the  onset  of  strain~rate  reversal  at  the  convex  outer  fiber 
(at  z  ■  a  when  die  >  0),  consider  the  strain  and  stress  Increment 
distributions  In  Fig.  3  for  die  >  0  associated  with  a  time  Increment 
dt.  Strain  rate  reversal  occurs  In  the  region  -a  <  z  <  z^,  so  elastic 
unloading  occurs  there.  From  Fig.  3a  showing  the  strain  Increments,  we 
have 


de  "  deQ  +  zdie  (2.7a) 

When  we  set  de  "  0,  equation  (2.1a)  determines  the  value  z^  that  divides 
the  cross  section  Into  loading  and  unloading  regions.  In  terms  of 
velocities 


z,  ■  -  u  /d» 

1  XX 


(2.8a) 


From  Fig.  3b  showing  the  stress  Increments,  we  have 

(H  de 

! E  de  -  a  <  z  <  z. 


do 


•i: 


*,<  z  <  a 


(2.9a) 
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By  Integrating  (2.9a)  with  de  given  by  (2.7a)  over  the  rod  eroaa 
section,  we  can  obtain  the  Incremental  thrust  and  bending  moment  in  the 
form 


**t  -  Hq  A  +  HlN/A  I  (2.10a) 

Mt  -  y/Tl  +  H2  I  o>x  (2.10b) 

where 

H  -  H  +  (E  -  H)i 
o  o 

Ht  -  (E  -  H)i1 

H2  -  H  +  (E  -  H)i2  (2.10c) 

*o  "  r  f‘l  6  * 

i.  bade 

vA  I  •'-a 

i  rzi  2 

i2  mf  J  b*  da  -  a  <  at  <  a  iet  >  0  (2.10d) 

when  *t  <  0,  the  thrust  and  moment  rates  are  still  given  by  (2.10a)  and 
(2.10b)  provided  we  change  the  lower  and  upper  integration  limits  in 
(2.10d)  to  aj  and  a.  For  a  circular  cross  section,  we  have  b  "  2(a^  “ 
s2^1/2  so  tj,e  integrals  in  (2.10d)  become 

i  -  -  ( f  .  -  4  sin  29, ) 
o  *1  2  1' 

*1  "  "  35  sln  h 
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(2.11a) 


i2  “I  (*i  ~78in  4*i>  Kt  >  0 

where 

cos  ■  -  z^/s  0  <  <  n  (2.11b) 

When  <  0,  we  replace  i0>  1^,  and  i2  by  j0>  J1#  and  j2  where 

jD  -  1  -  iQ  jl  "  “  il  J2  -  1  “  i2  (2.12a) 

In  Pig.  3a  (ic^  >  0),  we  see  that  partial  unloading  starts  when 
•  -a,  so  whenever  z  <  -a  the  Integrals  are  set  equal  to  aero.  In 
these  cases  (2.10c)  gives  HQ  -  H,  ■  0,  H2  -  fl  so  the  thrust  and 
moment  rates  (2.10a)  and  (2.10b)  reduce  to  (2.4a)  and  (2.4b)  for  the 
fully  loading  section.  If  partial  unloading  continues  with  2^  mono* 
tonlcally  increasing  until  it  reaches  the  value  2^  -  a,  we  reach  the 
stage  of  fully  elastic  unloading.  When  •  a,  we  have  i0  ■  1,  1^  •  0, 
and  i2  -  1  from  (2.10d)  and  consequently  HQ  «  E,  ■  0,  and  H2  -  E  from 
(2.10c)»  Then  the  thrust  and  moment  rates  reduce  to  (2.4a)  and  (2.4b) 
with  H  taking  the  special  constant  value  E. 

The  governing  equations  that  Include  monotonlc  partial  unloading 
are  the  set  (2.6a)  with  the  thrust  and  moment  rate  equations  replaced  by 
the  more  general  forms  (2.10a)  and  (2.10b).  In  matrix  form,  we  have 

A  w.  +  B  w  -  C  w  (2.13a) 

t  x 

In  which 
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0  -1 


-H0A  0  -HjAI  0 


-H]AI  0  -fl2I 


0  0 


0  0 


0  0 


0  0 


0  0 


0  0  0 


0  0  0 


0  0  1 


0  0  0 


0  0  0 


0  -k'AG  0 


0  0  -k’AG  000 


3.  CHARACTERISTIC  PROPERTIES  OF  THE  EQUATIONS 

Th  governing  equations  (2.13a)  form  a  system  of  quasi-linear 
hyperbolic  equations  of  first  order  (2.1b).  The  characteristic 
velocities  are  the  six  roots  of  the  characteristic  equation 


| A*1  B  -  c|  -  0 


Expanding  the  determinant  (3.1a)  leads  to 


(3.1a) 


[(pc2)  -  (Ho  +  H2)pc2  +  (HH2  -  H12)](pc2  -  k*G)  -  0  (3.2a) 


which  has  as  roots  the  six  wave  velocities 


cf  -± 


fof  k  +  H2  +  K  •  ®2>2  +  VI 


(3.2b) 


=.  ■  *  (h)'1 1".  +  "2  -  |(Ho  -  *2>2  +  -iff  (3-Jc> 


■  *  (# 


(3. 2d) 


The  fast  and  slow  wave  speeds,  and  c8  of  (3.2b)  and  (3.2c) 
govern  the  propagation  of  the  thrust  and  aoaent  In  both  directions  along 
the  rod.  The  wave  speed  c^  given  by  (3. 2d)  governs  the  propagation  of 
shear. 

Have  velocities  associated  with  a  fully  loaded  cross  section  are 
found  by  setting  H0  -  H2  ■  H  and  Hj  ■  0  in  (3.2b)  and  (3.2c).  These 
substitutions  give 

cf  "  cs  "  c  ”  <H/P>1/2  (3.3a) 

If  we  deteralne  the  left  null  row  vector  of  the  characteristic 
aatrlx  containing  the  corresponding  characteristic  value  c^,  that  la.  If 
we  solve 


(A_1B  •  Cj  I)  ■  0 


we  find  the  characteristic  equations  from 


*l(wt  +  C1  "  AiA_lc  w 


(3.4a) 


Following  this  procedure  leads  to  the  set  of  equations 


PCj  ARf  du  +  Rf  dM  +  pc#  I  d(i>  +  dM  -  ct  Q  dt 


£  *  WW  ■  SU  i  —  w  £ 

along  -  +  cf 


pc#  AR#  du  +  R#  dM  +  pc#  I  du  +  dM  -  c#  Q  dt 

along  ^  -  t  S 


pCq  A  dv  +  dQ  ■  +  k'AG  (1)  dt 
along  ^  -  +  c 


(3.5a) 
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where 


H  (I/A)17* 

Rf - * — T 

■*f  "Ho 


H^(I/A)'* 
pc  2  -  H 

8  O 


(3.5b) 


Characterlsclc  equations  associated  with  a  fully  loaded  cross 
section  are  found  by  setting  Hfl  ■  Hj  ■  H  and  Hj  -  0  In  the  governing 
equations  (2.13a)  finding  the  roots  (3. 2d)  and  (3.3a)  of  the  resultant 
determinant  (3.1a),  and  hence  the  resultant  null  left  row  vectors  of  the 
characteristic  matric  for  use  In  (3.4a).  Another  procedure  Is  to  find 
the  limits  of  (3.5b)  The  values  of  R^  and  Rg  tend  to  different  limits 
as  Cf  and  cg  tend  to  c.  Hence  the  equations  along  the  characteristics 
become 


pc  A  du  +  dN  ■  0 


pcIdu>  +  dM*cQdt 

dx 

along  gp  ■  +  c 

pc  Ad  +  dQ  -  +  k'  A  G  w  dt 
q  v 

along  -  +  cq 


(3.6a) 


4.  NUMERICAL  SCHEME 

To  Integrate  the  appropriate  form  of  the  characteristic  equation 
(3.5a)  along  the  characterlatlcs,  we  first  chose  a  rectangular  mesh  in 
the  x-t  plane  consisting  of  intersections  of  the  lines  x  ■  (1  -  l)Ax  and 
t  ■  (1  -  1) At  for  1  »  1,2...  and  J  ■  1,2....  He  then  adopted  a  scheme 
for  finding  the  values  of  the  dependent  variables  at  the  mesh  points  P 
along  the  line  t  “  tj^  from  values  Q  already  determined  along  line  t  • 
tj.  Figure  4  shows  typical  cells  with  linear  segments  of  the  backward 
and  forward  characteristics  required  to  determine  values  at  P  from  the 
values  of  Q_,  Q,  and  Q+  by  solving  the  finite  difference  form  of 


equations  (3.5a).  The  cell  size  was  determined  by  the  elastic  velocity, 

1/2 

that  is,  dx/dt  -  c  where  c  “  (E/p)  .  The  values  along  the  line  t  •  tj 

at  the  points  of  intersection  with  the  characteristics  from  P  were 
obtained  from  the  mesh  point  values  at  Q_,  Q,  and  0+  by  linear 
interpolation,  the  initial  slopes  c^  and  c8  of  the  plastic  wave 
characteristics  from  P  being  approximated  by  those  at  Q.  After  solving 
the  finite-difference  form  of  the  characteristic  equations  (3.5a)  to 
obtain  the  first  estimate  of  the  values  at  P,  the  dependent  variable 
values  at  the  midpoints  of  the  characteristic  segments  were  used  to 
adjust  the  slopes  for  a  second  iteractlon.  At  each  time  step,  z^  from 
(2.8a)  was  calculated  to  determine  the  local  loading  state  and  hence  the 
appropriate  form  of  the  characteristic  equation  could  be  selected.  If 
reloading  was  detected  by  a  reversal  of  the  monotlc  progress  of  Zj,  the 
calculations  were  terminated  because  the  theory  does  not  treat  this 
case. 


5.  NUMERICAL  EXAMPLES 

Calculations  were  performed  for  a  rod  having  a  radius  of  1  cm,  and 
a  length  of  20  cm.  He  used  the  simplified  stress-strain  relationship 
shown  by  the  curve  in  Figure  5.  The  curve  consists  of  a  linear  elastic 
part  smoothly  joined  to  a  linear  work  hardening  part,  and  it  is  based  on 
a  static  tension  test  on  uranium  containing  3/4%  titanium  [7].  The 
corresponding  value  of  the  material  properties  are: 


Young's  modulus 

E 

■1 

160  GPa 

shear  modulus 

G 

- 

59.3  GPa  (using  v 

hardening  modulus 

R 

- 

35  GPa 

yield  stress 

*7 

- 

730  MPa 

density 

P 

m 

18.1  g/cm3. 

He  used  the  bilinear  loading  history  of  Fig.  6  at  one  end  of  the  rod 
while  maintaining  the  other  end  free  of  loading.  Pure  thrust  loading 
gives  yielding  at  Ny  ■  2.3  x  103  N.  This  loading  history  was  chosen  in 
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a  simple  form,  suitable  for  evaluating  the  performance  of  the  calcula- 
tlonal  procedure,  rather  than  in  a  form  based  on  experimental  or  con¬ 
tinuum  mechanics  codes  results.  Also,  performance  evaluation  mas 
restricted  to  the  Initial  wave  transit. 

Before  generating  results  for  the  complete  rod  problem,  we  compared 
the  results  generated  by  our  computer  program  with  the  analytic  results 
for  the  case  of  pure  elastic-plastic  thrust  (normal  Incidence).  In  the 
computer  code,  the  bilinear  stress-strain  curve  between  points  1  and  2 
In  Fig.  5  was  replaced  by  a  smooth  curve.  In  the  analysis,  this  region 
was  replaced  by  a  straight  line.  For  this  problem,  a  ramp  loading  of 
0.786  x  10^  N/psec  was  applied.  Figure  7  shows  the  distributions  of 
thrust  along  the  rod  when  the  applied  thrust  has  reached  2.3  times  the 
yield  thrust,  N^.  For  excellent  agreement,  the  code  curve  should  pass 
through  point  2  as  well  as  point  1.  Because  of  the  different 
representations  of  the  stress-strain  curve  between  points  1  and  2  In 
Fig.  5,  the  solutions  are  different  between  these  points  In  Fig.  7.  The 
plastic  wave  from  the  loaded  end  should  pass  through  point  2.  However, 
propagation  appears  faster  than  It  should  because  the  Interpolation 
procedure  la  performed  out  within  an  x-t  mesh  slsed  for  the  elaetlc  wave 
propagation.  Figure  7  gives  a  qualitative  measure  of  the  error 
Involved. 

A  similar  comparison  was  made  for  propagation  of  bending  moments. 

A  ramp  shear  load  was  applied  at  one  end  of  the  rod  and  the  code,  with 
the  yield  stress  set  to  sero,  was  arranged  to  produce  plastic  moment 
propagation  within  an  x-t  mesh  slxed  for  elastic  wave  propagation  along 
dx/dt  "  (E/p)*^.  The  result  was  compared  with  elastic  wave  propagation 
by  changing  the  mesh  slse  to  dx/dt  “  (H/p)*^.  A  comparison  of  the 
moment  distributions  along  the  rod  Is  provided  by  Fig.  8.  The  effect  of 
the  Interpolation  for  plastic  waves  with  an  elastic  wave  mesh  la  seen 
with  the  solid  line.  It  Indicates  that  a  signal  is  being  allowed  to 
propagate  faster  than  Its  true  velocity,  as  In  the  case  of  plastic 
thrust  (Fig.  7). 
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Figure  9  shove  the  distributions  of  the  thrust  and  moment  along  the 
rod  when  the  elastic  wave  has  traversed  the  entire  length  of  the  rod 
(20  cm).  The  stress-strain  curve  Is  that  of  Fig.  5  and  the  thrust 
loading,  applied  at  an  angle  of  2*  to  the  rod  axis.  Is  that  of  Fig.  6. 

He  see  that  a  bending  moment  Is  developing  near  the  loaded  end,  so  under 
the  Ideal  loading  conditions  of  this  example  the  end  could  soon  break 
off. 


In  this  example,  the  code  operated  In  a  simple  mode  with  fully 
elastic-plastic  loading  throughout.  Similar  examples  at  larger  angles 
of  attack  (~8°)  cause  partial  unloading  of  the  rod  cross  section,  but 
for  the  few  examples  investigated  this  partial  unloading  was  soon 
followed  by  reloading,  for  whl<ch  the  code  mechanics  do  not  apply. 

Partial  unloading  Increases  the  section  bending  stiffness  and  hence  the 
resistance  to  the  bending  nomeat  Associated  with  this  Increased 
stiffness  Is  an  Increased  wave  velocity,  so  that  neighboring  sections 
become  influenced.  Further  work  Is  needed  to  define  the  limits  of 
applicability  of  the  program  In  terms  of  the  problem  parameters,  such  as 
angle  of  attack. 


6.  CONCLUSIONS 

Our  main  conclusion  Is  that  the  structural  mechanics  formulation 
works  well  for  analyzing  the  elastic-plastic  response  of  rods  under 
thrust  loadings  inclined  at  small  angles  to  the  rod  axis.  The  angles 
allowable  will  depend  on  other  parameters  such  as  rod  radius,  but  if  the 
bending  moment  Is  stimulated  as  a  growing  perturbation  on  the  thrust, 
useful  results  will  be  obtained  by  a  simple  computer  code. 

Further  numerical  Investigations  are  needed  to  extend  the  code  so 
that  useful  results  may  be  obtained  when  partial  unloading  occurs.  Our 
examples  Indicate  that  partial  unloading  Is  rapidly  followed  by  reload¬ 
ing.  Reloading  requires  monitoring  of  the  strain  history  and  Introduces 
a  complication  that  would  remove  the  advantage  of  simplicity  relative  to 
a  finite-difference  code. 
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(a)  Panetrator  under  almost  normal  impact 


(b)  Compressive  stress-strain  curve 


JA-1082-4 


FIGURE  1  ROD  CROSS  SECTION  STRESS-STRAIN  STATES  DURING  LOADING 


h*H 


(b)  Kinematics  of  a  rod  alamant 


JA-1 062-6 


FIGURE  2  GENERALIZED  FORCES  AND  DISPLACEMENTS  FOR  A  R00  ELEMENT 
(Sign  convention:  all  quantities  shown  positive) 
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(b)  Partial  unloading 


FIGURE  4 
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CHARACTERISTIC  CURVES  AND  MESH  POINTS 


TIME  T,  (js 

JA-1062-9 

FIGURE  6  BILINEAR  LOADING  FUNCTION 
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AXIAL  FORCE  RATIO.  -N/Ny 


DISTANCE  ALONG  ROD,  x  —  cm 

JA- 1062-10 

FIGURE  7  COMPARISON  OF  THRUST  RATIOS  FROM  ANALYSIS  AND  COMPUTER  PROGRAM 
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THRUST  AND  MOMENT  RAT 


FIGURE  9  THRUST  AND  MOMENT  DISTRIBUTIONS  AT  ELASTIC  WAVE  TRANSIT  TIME 
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